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Abstract 

This article is the continued version of the analytical blowup solutions for 2-dimensional 
Euler-Poisson equations in [10] and [11]. With the extension of the blowup solutions with radial 
symmetry for the isothermal Euler-Poisson equations in R^, other special blowup solutions in 
with non-radial symmetry are constructed by the separation method. 
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1 Introduction 

The evolution of a self-gravitating fluid (gaseous stars) can be formulated by the isentropic Euler- 
Poisson equations of the following form: 

{Pt+V ■ (pu) =0, 
(pM)t+V • (pw® u) + VP =-pV$, (1) 
A$(t,f) =a(iV)p, 

where a{N) is a constant related to the unit ball in : a{l) — 2; a{2) = 2tt and For > 3, 

aiN) ^ NiN - 2)ViN) ^ jV(jV ~ 2) ^^^^^ ^ , (2) 

where V{N) is the volume of the unit ball in and F is a Gamma function. And as usual, 
p — p{t,x) and u = u{t,x) = (ui, U2, Mat) 6 are the density, the velocity respectively. 
P = P{p) is the pressure. 
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In the above system, the self-gravitational potential field $ = ^{t,x) is determined by the 
density p through the Poisson equation. 

The equation (1)3 is the Poisson equation through which the gravitational potential is deter- 
mined by the density distribution of the density itself. Thus, we call the system (1) the Euler- 
Poisson equations. The equations can be viewed as a prefect gas model. The function P = P{p) 
is the pressure. The 7-law can be applied on the pressure P{p), i.e. 



Pip) = Kp'^:=!^, 

7 



(3) 



which is a commonly the hypothesis. The constant 7 = cp/c„ > 1, where cp, Cy are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in (3). In particular, the fluid is called isothermal if 
7 = 1. It can be used for constructing models with non-degenerate isothermal cores, which have a 
role in connection with the so-called Schonberg-Chandrasekhar limit [7] . And we denote the radial 
diameter as: r := 



The system can be rewritten as 



Pt 



(4) 



V ■up + Vp-u={), 
P + Er=i ^k^) + £-P{p) —p£-Hp), for ^=l,2, ...N, 

A<^{t,x) = a{N)p. 

For N = 3, (4) is a classical (non-relativistic) description of a galaxy, in astrophysics. See [2], [3] 
and [7] for a detail about the system. 

For the local existence results about the system were shown in [1] and [5]. Historically in 
astrophysics, Goldreich and Weber constructed the analytical blowup (collapsing) solutions of the 
3-dimensional Euler-Poisson equations for 7 = 4/3 for the non-rotating gas spheres [6]. After that, 
Makino [8] obtained the rigorously mathematical proof of the existence of such kind of blowup 
solutions. Besides, Deng, Xiang and Yang extended the above blowup solutions in R'^ {N > 4) 
[4]. Recently, Yuen obtained the blowup solutions in R"^ with 7 = 1 by a new transformation [10]. 
The family of the analytical solutions with radial symmetry 



u=-V{t,r), 



are rewritten as 

For TV > 3 and 7 = (27V - 2)/N, in [9] 
' 1 

p{t, r) = < 



-,N 



for r < a{t)Z^; 



N - 1 



a{t)-- 



0, for a{t)Zf, < r. 

a(0) ^ 0, a(0)= 02 



(5) 



{2N - 2)K 



N/{N-2)_ 



= M, y(0) = a > 0, y(0) = 0, 
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where /x = [N{N - 2)X]/{2N - 2)K and the finite is the first zero of y{z); 
For = 2 and 7 = 1, in [10] 



p{t,r) = 



1 



oV{r/a(t)) 



a{t)=—T-, a(0) = > 0, d(0)= 02; 

a[t} 

a{2) y^'^ 



y{z)+-Jj{z)+ ^ 



II, 2/(0) = a, 2/(0) = 0, 



(6) 



where K > 0, jj. = 2X/K with a suflaciently small A and a are constants. 

However, all the known solutions are in radial symmetry. In this paper, we are able to obtain the 
similar results to the non-radial symmetric cases for the 2-dimensional Euler-Poisson equations (4) 
in the following theorem. 



Theorem 1 For the isothermal Euler-Poisson equations (4) in , there exists a family of solu- 
tions, 



„2_| 2 N 

a{t) 

— +'^, u{t,x) 



a{t) 



2a{t) 

a{t) = ^ia2t + 2aif 



{xi,X2,xi +a;2,...a;i -|-a;2). 



<I.(s) 



s$(s) + ^ - e*e" K = 0, $(0) = a, $(0) = e*e~K, 



(7) 



where A, B, ai 7^ 0, 02 , 



a{N)e'' 



e* > 0, a and (3 are constants. 



In particular, oi > and 02 < 0, the solutions (7) blow up in the finite time T = —02/01. 



2 Special Blowup Solutions I 



Before presenting the proof of Theorem 1, we prepare the following two lemmas. 



Lemma 2 For the continuity equation (4)1 in R , there exist solutions. 



p{t, x) 



a{t) 



, u{t,x) 



a(t) 
2a{t) 



{xi,X2,X-i_ -\-X2,...,Xi -\-X2), 



(8) 



where the scalar function f{s) > € and a{t) 7^ G C^. 
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Proof. We plug the solutions (7) into the continuity equation (4)i, 



pt + 'V ■ up + Vp ■ u 



d_ 

dt 



+ v- 



a{t) 



+ V 



2a{t) 



f 



{xi,X2,Xi +X2,:;Xi +X2)- 



a{t) 



a{t) 



a(t) I ait) , , 
- -77r(a;i, xi, xs + X2, ....a;2 + a;2) 



a{t) a{t) 



r f ^1 + 1 d ( x\ + X 



(9) 
(10) 

(11) 



a{t) I d 



d 



N 



d 



^2 + ^2 

a(i) 



«(*) 



2a{t) 



d 



f 



dxi a{t) 



■ Xl + 



d 



f 



8x2 a{t) 



^ d ^ 



i=3 



'dXi 



a{t) 



+ 



a?{t) ' V J ait)' \ a{i) 

d{t) 



I 



1 >W^V «W 



,,.2 



a{t) a?{t) 



2a{t) 



. xl + X2 \ i [ x\~^ ^2 



a{t) j 2x1 V / 



= 0. 



(12) 



[xi + a;2) 



(13) 

(14) 

(15) 
(16) 



The proof is completed. ■ 

The following lemma handles the Poisson equation (4)3 for our solutions (7): 



Lemma 3 The solutions, 



1 "\~sw , 

o=—e- ^ +^ (17) 



with the second-order ordinary differential equation: 



^(s) Mf) • a 

s$(s) + - e*e-~ = 0, $(0) = a, $(0) = e*e-K, (18) 



w/iere s := {xf + X2)/a{t) and C, "^^J^ = e* and a are constants, 
fit into the Poisson equation (4)3 in R'^ . 
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Proof. We check that our potential function ^{t,xiX2) satisfies the Poisson equation (4)3: 



A^{t,x) - a{N)p 



-l+xl\ a{N) 



5W 



+c 



{t) J ait) 



= V- 

_ a{N) 

= V- 

- A. 

dxi 



dxi \ a.{t) J ' dx2 \ a{t) J ' dxs \ a{t) 



d ^( xl+xl 
a{t) 



dx 



JV 



5Ttr 



+c 



^ xl+xl \ ^ xl+xl \ 2X2 Q 



a{t) J a{t) V a{t) J a{t) 
d 



+ x\ 2xi 
^ a{t) 'a{t) 



+ 



dx2 



■ xj+xl 2x2 
ait) 'a{t) 



ait) 

aiN) 



^ o(t) ^a2(t) ^ ^ ait) 'ait) 

2 I ^2 ^^2 ^2 I ^2 o ^JAT\ * -RtT 



^ xj+xj 4xi - xj+x^ 2 aiN) 
^ ait) U^it) ^ ait) 'ait) ait) 

4 ( ■■ , , $(s) aiN)e^ 



where we choose s := ix\ + x5)/a(t) and the ordinary differential equation: 
s$(s) + - e*e-~ = 0, $(0) = a, $(0) = /?, 



(19) 

(20) 
(21) 

(22) 

(23) 

(24) 
(25) 

(26) 
(27) 



(28) 



with '^^^^'^ = e*^ a and /? are constants. Therefore, our solutions (71) satisfy the Poisson equation 
(4)3. 

The proof is completed. ■ 

Besides, we need the lemma for stating the property of the function $(s) of the analytical 
solutions (18). We need the lemma for stating the property of the function $(s). In particular, 
the solutions (7) in A/'-dimensional case involve the following lemma. The similar lemma was 
already given in Lemmas 9 and 10, in [10], by the fixed point theorem. For the completeness of 
understanding the whole article, the proof is also presented here. 



Lemma 4 There exists a sufficiently small e* > 0, such that the ordinary differential equation 

— n 

(29) 



s^s)+^-s*e K =0, 



$(0) = a, $(0) = -£*e K 
where K > and a are constants, has a unique solution $(s) G C^[0, 00). 
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Proof. The lemma can be proved by the fixed point theorem. The equation (29) can be rewritten 
as: 



= 1/2. 



ds 



ds 



(sl/2$(5) 



(sV2$(a)) = 



_*(£) 

-s*e K , 

*(£) 

-£*e" K 



,1/2 



(30) 
(31) 



With the initial conditions: $(0) = a and $(0) = — e*e k , the equation (29) is reduced to 

4>(t) 

$(a) = / ^rfr. 



el/2 / T-1/2 
S ' Jo ^ ' 



(32) 



Set 



a.(r) 



For any so > 0, we get / G C^[0, so]- And for any <l>i, $2 S ^^[0, sq], we have, 



(33) 



|/(s,<i>i(s))-/(s,$2(s))| 



s I e~ K -e~ K 



dr 



,1/2 



(34) 



As e* is a function of we can show that the function e*, is Lipschitz-continuous. Then we 
get, 



|/(s,$i(s))-/(s,$2(s)| = 



(<I>2(r)-$i(T) 

fT72 



dT 



<- sup |<I>i(T)-$2(r)|. (35) 

0<s<so 



Let 



r$(s) = a+ [ J{T,^{T))dT. 
Jo 



(36) 



We have T$ e C[0, sq] and 



|T$i(s)-T$2(s) 



r .f{T,^^{T))dT - r f{T,^2{r))dl 
Jo Jo 



<- sup |$(s)i-$(s)2|. (37) 

0<s<so 



By choosing the constant e* suck that < ^ < 1, this shows that the mapping T : C[0, sq] — > 
C[0, Sq], is a contraction with the sup- norm. By the fixed point theorem, there exists a unique 
$(s) e C[0, So], such that T$(s) = $(s). 

It is because that the chosen constant e* is independent of the variable x. Therefore, we have the 
global unique solution $(s) G C^[0, 00). The proof is completed. ■ 

Now, we are ready to check that the solutions fit into the Euler-Poisson equations (4). 
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Proof of Theorem 1. By Lemma 7 and Lemma 8, the solutions (7) satisfy (4)i and (4)3. For 
the Xj-component of the isothermal momentum equations (4)3 in (-/V > 3), we have 



N 



f dui ^ dui j _|_ i9 ^ _|_ 5^ 



fc=l 



dx 



= P 



+ K 



d a{t) I d{t)xi g a(t)xi a(t)x2 Q a(t)xi 

dt 2a(t) -^l ' 2a{t) dxi 2a (t) 2a(t) dx2 2a(t) 

I a{l)(xi+X2) d ( a(t)xi \ 

2a(t) dxi \ 2a(t) J 



d e- 



i=3 

1 



+c 



dx a{t) 

1 /a(i) _ 

2 \a{t) a?{t) 



X-i ~\~ Xn 



a{t) J a{t) 



1 a{t)xi d{t) 

•^^ 1 ~l~ "T 



4 a{t) a{t) 



2,2 



Jit) 



- K- 



a(t) J Ka{t) ^ 



= P 



= 0, 



1 /2'd{t) r?{t) 
4 \ ' 



Xl 



a(i) y a{t) 



a{t) 
2xi 



2x1 



a{t) 



2xi 
W) 



(38) 
(39) 

(40) 
(41) 

(42) 

(43) 
(44) 



where we used 



2a{t)a{t) - a\t) = 0, a(0) = ai > 0, d(0) = 02- 



(45) 



which is exactfy solvable by Maple, 



a{t) = ^(o2^ + 2al)^ 



(46) 



For the a;2-component of the isothermal momentum equations (4)3 in R" , we have 



/ du2 , du2\ 



+ —Kp + p— 
0x2 0x2 



d a(t)x2 d{t)xi d d{t) d{t) d d{t)x2 



N 



dt 2a{t) ^ 2a{t) dxi 2a{tf^ ^ 2a{tf^ 8x2 2a{t) ^ '^'^^ Oxu 2a{t)^^^ ^ 



d d{t) 



+ K 
= P 

= P 
= 0. 



d e 



dx a2(i) 
1 fd{t) a^{t) 



hp 



2 \a{t) a?{t) 
1 f2d{t) a^it) 
4 \ ~ 



X2 + 



X2 



dx2 
d{t)x2 d{t) 
Aa{t) 'aif) 



+ + 



(47) 
(48) 

(49) 

(50) 

(51) 
(52) 
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For the a;j-component {i > 3) of the isothermal momentum equations (4)3 in R" , we have 



CdUi ^ ■y^ dUi j _|_ ^ ^ _|_ 9^ 
dt ^ dxk J dxi dxi 



d d(t)(a;i+a;2) i a(t)xi d a(t)xi . a(t)x2 d a(t)x2 
dt 2a{t) 2o(t) 9x1 2a(t) ~'~ 2a(t) dx2 2o(t) 

,,, d a(t){xi+X2) 
"I" 2^k=3 "fc dxk 2a(t) 



a(t) ^ 

+ — 



= P 

= P 
= 0. 



1 fa{t) a^{t) 



2 \a{t) a?{t) 
1 /2a(t) d2(t) 
4 ~ 



dxi \ a{t) 

a{t)xi d{t) d{t)x2 d{t) 

[Xl + X2) -\ — 7-^ — — 777 + 



2a{t) 2a{t) 2a{t) 2a{t) 



{xi + X2) 



(53) 
(54) 

(55) 
(56) 
(57) 
(58) 



Therefore, our solutions satisfy the Euler-Poisson equations. In particular, ai > and 02 < 0, the 
solutions (7) blow up in the finite time T = —02/01. 
The proof is completed. ■ 

It is clear to see the blowup rate of the solutions (7): 
Corollary 5 The blowup rate of the solutions (7) is, 



l\mp{t,6){T-ty>0{l). 



(59) 



3 Special Blowup Solutions II 

In the recent paper [11], we have the special solutions for the isothermal Euler-Poisson equations 
in R^, in the following from; 



I \ <»(*) / 



p{t,Xi,X2) = K +'^, u{t,Xj^,X2) = ^{X1,X2), 



a{t) = ai + a2t, 



(60) 



$(s) - e*e-^ = 0, $(0) = a, 4>(0) = /?, 



where A, B, ai 7^ 0, 02 , ^2^g2 = e* > 0, a and /? are constants. 

In this section, we extend the above solutions to the A/'-dimensional Euler-Poisson equations 
(4) in the following theorem: 



Theorem 6 For the isothermal Euler-Poisson equations (4) in R {N >3), there exists a family 
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of solutions, 



'°(*'^) = ^ ^{t^x) = —{xi,X2,Xi,...,Xi) 



ait) 



a{t) = ai+ a2t, 



(61) 



$(s) - e*e~K^ = 0, $(0) = a, $(0) = /3, 

where A, B, ai 7^ 0, 02 , = e* > 0, a aKC? /3 are eonstants. 

In particular, ai > a^d 02 < 0, i/ie solutions (61) blow up in the finite time T = —02/01. 



Before presenting the proof of Theorem 1, we prepare the following two lemmas. 
Lemma 7 For the continuity equation (4)1 in , there exist solutions. 



f 



p{t,x) 



ait) ) «(0 / \ 

-^ij^ , U(t,x) = ^^{Xi,X2,Xi,...,Xi), 



(62) 



where the scalar function f{s) > e and a{t) 7^ G C^. 

Proof. We plug the solutions (61) into the continuity equation (4)i, 

Pt + V • up + V/9 • M 



d_ 

di 



+ v- 



I Axi+Bx2 
^ \ <t) 
aHt) 

Axi + Bx2 



_ d{t) , 
+ V • — — (a;i,a;2,a;i, ....,x\ 
a{t) 



d(t) 

— {xi,X2,Xi,...,Xi) 

a[t) 



f 



Ax\ + Bx2 



~2a{t) ^ f Axi + Bx2 



a{t) 



a^(t) dV 



N 



1 d jAxi+Bx2\ hit) d d 

' a{t) ) + W) \d^^"' ^ 5^"' ^ h'^^" 



(63) 
(64) 

(65) 



Axi + Bx2 

ait) 
a^t) 



(66) 



a{t) 



/ Axi + Bx2 \ 



dx\ 



a^{t) 



Xi + 



Ax\ + Bx2 



dxo 



'X2+J2 



N g f 



Axi + Bx2 



i=3 



dxi 



a^it) 



■Xi 



(67) 



-2a{t) ^ Axi+Bx2 
a^{t) ^ V a{t) 



1 ; ^ Axi+Bx2 \ (Axi + Bx2)d{t) _^ ^a{t) 



a2(t)- 



a{t) 



f 



Axi + Bx2 



a{t) a^{t) 



+ 



d(t) 
a{t) 



Axi + Bx2 



a2(t) a{t) 



■ / Axi + Bx2 
^xi ait) 



BX2 



a^it) 



ait) 



(68) 
(69) 



= 0. 



(70) 
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The proof is completed. ■ 

The following lemma handles the Poisson equation (4)3 for our solutions (61): 



Lemma 8 The solutions, 



n a(t) ) 



+c 



with the second-order ordinary differential equation: 



$(s) - e*e~~ = 0, $(0) = a, $(0) = /3, 



(71) 



(72) 



where s {Axi + Bx2)/a{t) and C , "^a^gs = e*, a and (3 are constants, 
fit into the Poisson equation (4)3 in . 

Proof. We check that our potential function $(t,xi,X2) satisfies the Poisson equation (4)3: 



A<i>{t,x^,x2)-aiN)p 

'Axi+Bx2\ a{N) ^_<z!^^^c 



V- 



a{t) J oH^) 



d ^f Ax.+Bx-A d f Ax,+Bx2 \ d ^/ Ax,+Bx2 
dxi \ a{t) J ' dx2 V o,{t) J ' dxs \ a{t) 



dx 



(73) 



(74) 

d ^fAj-.+Bx^ 



N 



a{N) K^^: 



V- 

d 



■ / Axi + Bx2 \ A ■ / Axi +Bx2 \ B ^ 
ait) J a{t) ' ^ ait) ) ait) ' 



.(TV) 



a^it) 



+c 



■ Axi+Bx2 . A 
a{t) a{t) 



+ 



d 

dX2 



■,Axi+Bx2. B 
a{t) 'a{t) 



,iN) _<=m=l 



+c 



A' + B^ 



$(s) 



aiN)e^ 

— — e K 

A^ + B^ 



where we choose s := {Axi + 8x2) /a{t) and the ordinary differential equation: 

$(s) - e*e-~ = 0, $(0) = a, $(0) = /3, 



a{t) 
(75) 

(76) 
(77) 
(78) 
(79) 



(80) 



with '^a^ga = e* , a and (3 are constants. Therefore, our solutions (71) satisfy the Poisson equation 
(4)3. 

The proof is completed. ■ 

Now, we are ready to check that the solutions fit into the Euler-Poisson equations (4). 
Proof of Theorem 1. By Lemma 7 and Lemma 8, the solutions (61) satisfy (4)i and (4)3. For 
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the Xi-component of the isothermal momentum equations (4)3 in {N >3), we have 



dui 
~dt 



N 

E 

k=l 



Uk 



dui 
dxk 



dxi dxi 



N 



Uk- 



d a(t)xi 



d a{t)xi a{t)xi d a{t)xi a{t)x2 d a{t)xi 
dt a{t) a{t) dxi a{t) a{t) 8x2 a{t) ' ^^'"'^ dxk a{t) 



+ K 



d e 



+c 



dxi 

d{t) a?{t) 



/ Axi + Bx2 \ A 
V <t) ) <t) 



a{t) a^{t) 



Xl + 



~K 



K- 

■ P 
■0, 



am 

d{t)xi 



d{t)xi d{t) 
a{t) 'm 

Axi + Bx2 



a{t) 

Axi + Bx2 \ A 



A ■ / Axi + Bx2 

Axi + Bx2 \ A 



A 

W) 



by taking d{t) =0 that is 



a{t) = ai + 02*- 



For the a;2-component of the isothermal momentum equations (4)3 in , we have 



N 



— + uk^^ I + —Kp + p— 

dt ^ dxk j 8x2 9x2 



+ K 



= P 



d d{t)x2 d{t)xi d d{t)x2 d{t)x2 d d{t)x2 

a{t) 8x2 a{t) 

) 



dt a{t) a{t) 8x1 a{t) 
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8x2 o?-{t) 
d{t) a^{t) 



X2 + 



8X2 

d{t)x2 d{t) 



N 

E 

fe=3 



Uk 



8 d{t)x2 
8xk a{t) 



-K- 



a{t) a'^{t) J a{t) a{t) 

■f Axi+Bx2 \ B .f Axi+Bx2 \ B 

V a{t) J Ka{t) ^ V <t) ) ait) 



am 

d{t)x2 
a{t) 
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For the Xi{i > 3) -component of the isothermal momentum equations (4)j in , we have 

N 



du. 



du3 



dt ' ^^^'^Sxfc j dx. 



k=l 



d d{t)xi a(t)xi d a{t)xi a{t)x2 d a{t)xi ■^-^ d a{t)xi 
dt a{t) a{t) dxi a{t) a{t) 8x2 a{t) dxk a(t) 

a(t) d'^{t)\ d{t)xid{t) 



a{t) a{t) 



a{t) 
a{t) 

= 0. 



(96) 

+ + (97) 

(98) 
(99) 
(100) 



Therefore, our solutions satisfy the Euler-Poisson equations. In particular, ai > and 02 < 0, the 
solutions (61) blow up in the finite time T = —02/01. 
The proof is completed. ■ 

Remark 9 The existence and uniqueness of the function ^(s) in the solutions (61)3 can be shown 
by the fixed point theorem, if e* is a sufficient small number. 

Additionally the blowup rate about the solutions is immediately followed: 

Corollary 10 The blowup rate of the solutions (61) is, 



\imp{t,0){T -ty > 0(1). 



(101) 



Remark 11 The other analytical solutions in can be constructed by 



K +'-,u{t,X) = ^-r{Xi,X2,X3,...,XN), 

a(t) 

a{t) = ai+ a2t, 



$(s) - e*e~R^ = 0, $(0) = a, $(0) = 



where Xi = xi or X2- 



(102) 



Remark 12 Our solutions (7), (61) and (102) also work for the isothermal Navier-Stokes-Poisson 
equations in R^ {N > 3): 



Pt+V ■ ipu) =0, 
{pu)^+V ■ {pu (Siu) + \/Kp =-/3V$ + fiAu, 



(103) 



A<^{t,x) = a{N)p, 



where fx > is a positive constant. 



In conclusion, due to the novel solutions obtained by the separation method, the author conjec- 
tures there exists other analytical solution in non-radial symmetry. Further works will be continued 
for seeking more particular solutions to understand the nature of the Euler-Poisson equations (4). 
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